NASA 

Technical 

Paper 

2252 


January 1984 




NASA fv 

TP I a 

2252 f:>. 




Data 



LOAN COPY: RETURN TO 
AFWL TECHNICAL LIBRARY 
KIRTLAND AFB; N.M. 87117 


Larry D. Lee, /■%’: 
George B. Finelli, 
Mitchel E. Thomas, 
and Felix L. Pitts 


_<- G 


_<p rP ^ <y g> 


J&\rSo O O ^ >%?*■ "' nJS’ a. 


W - ' -O. ■ • . -O . ^ *\S» .<C* . f? 


<A e,® 


<cr <>y*v ^ «6.- 


^>;*v 

^ <0-^ i* 


<fX<f <f 




NASA 

Technical 

Paper 

2252 

1984 


Statistical Analysis 
of Direct- Strike 
Lightning Data 
(1980 to 1982) 


Larry D. Lee, 
George B. Finelli, 
Mitchel E. Thomas, 
and Felix L. Pitts 

Langley Research Center 
Hampton , Virginia 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Office 

1984 


TECH LIBRARY KAFB, NM 

IlDHIHIII 

□ObflDQS 


INTRODUCTION 


The objective of the NASA Langley Research Center direct-strike lightning 
research is to aid in characterization of the direct-strike lightning environment 
affecting aircraft. The research is motivated by the need for refined characteriza- 
tion of the lightning-aircraft environment in support of advanced aircraft designs 
that contain composite structures and digital electronic systems. 

A specially instrumented F-106B aircraft is used in the research to penetrate 
thunderstorms and elicit direct strikes. Measurements of the lightning-aircraft 
electromagnetic interaction process are made with sensors derived from those designs 
developed for nuclear electromagnetic pulse measurements. Nine sensors with band- 
widths greater than 100 MHz are located about the aircraft and measure the rates of 
change of magnetic and electric flux density and current in the nose boom. Two 
transient recorders with bandwidths of 50 MHz each (100 MHz sample rates) record the 
output of two selected sensors, and the trigger configuration allows time correlation 
of simultaneously recorded waveforms. Two other transient recorders, triggered inde- 
pendently and operating at a sample rate of 60 MHz, record nose-boom and vertical 
tail-fin currents. An airborne field mill system gives a continuous indication of 
charge polarity on the aircraft as well as vertical and lateral electric fields from 
dc to 10 Hz. (A list of symbols used in this paper appears after the references.) 
Photographic time histories of lightning attachments to one wing and to the vertical 
tail have been acquired with a 1 6-mm movie camera with a wide-angle lens. A time- 
code generator enables time correlation of 1 msec between the photographic data and 
other recorded aircraft and ground-based measurements. Ground-based measurements 
include fast and slow electric field changes and radio frequency spherics generated 
by lightning in storms penetrated by the aircraft. 

The instrumentation system, transient recorders, sensors, and a summary overview 
of the tests are presented in references 1 to 5. Interpretation, analysis, and gen- 
eralization of the direct-strike data involves a methodology for deriving the light- 
ning source from the response measurements using electromagnetic coupling codes as 
described in reference 6 along with time-domain ref lectometry techniques using lab- 
oratory aircraft scale models for resonant frequency identification and comparison 
with flight data as described in references 7 and 8. The complete set of waveforms 
recorded in 1980, 1981, and 1982 are reported in references 9, 10, and 11. 

Data from 10 strikes to the aircraft were recorded during each of the first 
2 years of tests in 1980 and 1981. In 1982, data were recorded from 156 strikes. 

The significant increase in the number of strikes in 1982 was brought about by 
improved operational procedures. An aircraft fuel system modification provided an 
additional 1 /2-hour mission time for a total mission time of 1 1/2 hours, the pene- 
trations were made at generally higher altitudes of 25 000 to 35 000 ft, and guidance 
to specific altitudes of lightning activity was obtained based on UHF radar returns 
from lightning channel ionization. 

The measurements made with the direct-strike data system are nose-boom current 
I and its rate of change I, rate of change of electric and magnetic flux density 
D and B, and static electric field E. The nominal sensor locations are shown in 
figure 1. # The I and I sensors are located at the base of the nose boom at posi- 
tion 1; D sensors are located under the fuselage at position 2, under the wings at 



positions 9 and 11, and on the port side of the vertical fin at position 12 corre- 
sponding to the respective sign conventions for D R , D w L , R , and D^; B sen- 

sors are located at opposite sides of the aft fuselage at positions 6 and 7 and under 
the wings at positions 8 and 10 corresponding to the respective sign conventions 
for B l (starboard side), B T (port side), , and R ; and E sensors are 

located at position 3 corresponding to sign conventions for' E z and at positions 4 
and 5 corresponding to Ey. The sign conventions established for the measurements 
shown in figure 1 indicate positive quantities in the direction of the arrows? the 
relation between the sign conventions of the variables D, B, and I and their 
derivatives are also shown in figure 1 in the measurement polarity time history. All 
the data reported and analyzed in this paper are from the B sensor at position 6, 
the D sensor at position 2, and I and I sensors at position 1. 

The statistical analyses and results presented in this report provide insight 
into the expected ranges of the electromagnetic environment when fighter-class air- 
craft are struck in thunderstorms at altitudes of 20 000 to 35 000 ft. It should be 
noted that the statistics characterize the first transient of the lightning process 
which is greater than the system threshold and which occurs within the first 
1300 usee of the time the threshold is exceeded. Data acquired subsequently but too 
late for inclusion in the 1982 statistical analysis data base indicate these statis- 
tics are actually representative of the entire lightning process. It is known from 
limited bandwidth continuous recordings that the total lightning process consists of 
a moderate number of discrete electromagnetic pulses distributed over a second or so 
(ref. 7). The continuous records indicate the transient recorder threshold is suf- 
ficiently high so that the time of its exceedance is distributed over the duration of 
the lightning process (and not limited to particular time intervals in the process), 
and wide-band peak detectors which monitor the peak amplitudes of the entire process 
yielded results commensurate with the 1982 data statistical predictions. 

The primary objectives of the statistical analysis are: to identify a location- 

scale family of distributions consistent with the data, and to use the resulting 
model to compute estimates of the quantiles. Graphical display of the data suggests 
a lognormal distribution for some subsets. Statistical significance tests support 
the hypothesis that the D and I measurements are lognormally distributed. For 
these data, and for the B measurements, quantiles and confidence limits are 
established. 

The cumulative direct-strike lightning data base includes all the waveforms 
recorded during 3 years of flight tests (1980 to 1982) with transient recorders. The 
recorded data are categorized as resulting from direct lightning strikes or from 
nearby lightning events. If the lightning attachment to the aircraft was observed by 
the flight crew or by the onboard cameras, or if the current sensors indicated cur- 
rent flow, the lightning event was identified as a "strike." Other lightning occur- 
rences, which were not confirmed as strikes by the above criteria, were termed 
"nearby" events. The notation to identify data in these two categories for the sta- 
tistical analyses is as follows: B^ and identify B and D data for nearby 

. . • • . 

events, and B g and D g identify B and D data for strikes; I and I data 

result only from strikes. 
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STATISTICAL RESULTS 


The span (range) of the recorded waveforms was arbitrarily chosen as the random 
variable in the statistical analyses. The span is the difference between the maximum 
positive and the maximum negative amplitude of the measured waveforms. Table 1 shows 
the number of records acquired for each of the six distinct data sets (referred to as 
samples) and shows the range of span values for each sample. 

For each of the samples , frequency-of-occurrence tables were generated to sug- 
gest the underlying distribution of the random variable. (See tables 2 to 7 . ) For 
each of these tables, the range of the particular data set is broken into 10 equally 
sized intervals. The lower limit of each interval is shown in the first column. The 
second column shows the frequency of data values which are greater than or equal to 
that lower limit but less than the lower limit of the next interval. Samples 
and D n were excluded from further analysis since the variable span^is dependent 
on the distance to the lightning source, which was not known. The I sample was 
also excluded from the following analysis on the basis of insufficient sample size. 


Probability Plots 

The probability plotting method described in the appendix and known as the Q-Q 
plot (ref. 1 2) is used to judge the adequacy of the lognormal and the type II extreme 
value distribution of largest values. Since this plotting method is useful only for 
judging the adequacy of location and scale parameter models, the original sample 
values x 1 , x , . .., x^ are replaced by their natural logarithms, so that 

y^ = ln(x^) for i = 1, 2, . .., n. This transformation also converts the lognormal 
and type II extreme value distributions to the normal and type I extreme value dis- 
tributions, which are location and scale parameter models. 

Members of the set y^ , y^, ..., y n are ordered from smallest to largest, and 

the ith ordered value is paired with a measure of location of the ith ordered value 
in a sample from the theoretical distribution. In the plots, departures from linear- 
ity suggest that the lognormal and type II extreme value distributions are not ade- 
quate models. Tables 8 to 1 0 are also helpful in judging the adequacy of these 
distributions when compared with the theoretical distributions. 

Figures 2 to 4 show plots of the ordered transformed samples # versus the theo- 
retical quantiles of the standard normal distribution. For the D g and I samples 
(figs. 3 and 4), the plots are nearly linear, suggesting a good fit to the normal 
distribution. Thus, the lognormal distribution is an adequate model for the true 
parent distribution giving rise to the untransformed sample x^ , x 2 , ..., x n# In 

the case of the B g sample (fig. 2), the plot exhibits a concave shape and suggests 
that negative skewness is a property possessed by the true parent distribution of 
y 1 , y 2 , ..., y n . The significance test for normality, which is presented later, 

supports the conclusion that the concave shape is a significant departure from the 
linear shape expected under normality. 

Figures 5 to 7 show plots of the ordered transformed samples versus the theo- 
retical quantiles of the type I extreme value distribution of largest values. Gen- 
erally, none of these plots indicate that the fit to the extreme value model is 
better than the fit to the normal distribution model. However, the heavier upper 
tail (positive skewness) of the extreme value model and a possible link of this model 
with the physical feature that span values are themselves largest values motivated 
the a priori selection of the extreme value model. 
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A significance test for adequacy of the lognormal model can be based on the 
probability plot correlation coefficient (ref. 'll).' The small observed significance 
level (less than 0.005) shown in table 11 for the B g sample suggests a significant 
departure from the lognormal distribution model. As anticipated from the probability 
plots, the D q and I samples do not yield significant values of the correlation 
statistic, indicating the data support the hypothesis that the underlying distribu- 
tions are lognormal. 


Quan ti le Es tima tes 


Asymptotically unbiased quantile estimators as developed in reference 14 and 
described in the appendix are applied to estimate the upper quantiles y^ of the 
true parent distribution for y-j , y 2 , . .., y n * Then the transformation 

x^ = exp(y ) gives an estimate of the quantile x^ of the true parent distribution 
for x-j , x 2 , • Tables 12, 13, and 14 show the estimator x~ of the quantile 

x r , its es tima ted # variance, and the 95-percent lower and upper confidence limits L 


iffd U for the 


B c 


D 


g/ and I samples* 


Sample Size 

Of particular importance in direct-strike lightning experiments is the problem 
of estimating the "worst case" values of the distribution of the measurements. The 
extreme quantiles of most interest require large sample sizes to obtain "good" pre- 
cision for the estimates. Precision is commonly measured by the length of its asso- 
ciated confidence interval. However, because the theoretical confidence limits have 
an expected length depending on the unknown quantile, precision as measured by length 
would also depend on the unknown quantile. To avoid this problem, we examine the 
effect of sample size on the ratio R = U/L. A general representation of the sample 
size needed to ensure that (1 - a) 1 00-percent confidence limits have a ratio not 
exceeding a specified R q is given in the appendix. 

Table 15 shows how the sample size varies with the required precision as mea- 
sured by R. The computations were made for 95-percent confidence limits for x^ 
assuming that estimation is according to the normal distribution. For D q samples 
and for estimating x 95 / a sample of size n = 46 yields an R q = 2.0, whereas for 
estimating x ^ it takes a sample of size n = 1 44 to yield the same precision. 
Similar conclusions are obtained for other samples by examining table 15. 


CONCLUDING REMARKS 

Statistical methods have been applied in the analysis and interpretation of 
electromagnetic data acquired during 176 direct strikes to the F-106B aircraft. 
Probability plots and significance tests suggest that the data for the rate of change 
of electric flux density D and for the current I are lognormally distributed. 

The data for the rate of change of magnetic flux density B, on the other hand, 
appear more uniformly distributed. An explanation is that truncation occurred in the 
measurement process - some of the data points in the upper interval are at the full- 
scale limits of the measurement system, for example. 

# Techniques were presented for estimating the quantiles of the distributions of 
B, D, and I measurements. Hie estimate of the theoretical 90th quantile of the 
distribution yielding the B sample for direct lightning strikes is 1497.09 T/sec, 
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and 95-percent confidence limits e for this quantile range from 1075.83 T/sec to 
2083.31 T/sec. Similarly/ for D the estimate is 16.10 A/m 2 / with limits of 
13.23 A/m 2 and 19.58 A/m , and for I the estimate is 9156.66 A, with limits of 
7032.07 A and 11 923.15 A. The largest observed values for direct strikes were 
1727 T/sec, 30.60 A/m 2 , and 13 900 A for B, D, and I. An insufficient number of 
samples for rate of change of current I measurements for statistical analysis have 
been recorded to date. The maximum zero-to-peak I recorded, however, exceeded the 
full-scale peak of 48 GA/sec; the corresponding maximum I derived graphically from 
the leading edge of the simultaneously recorded current waveform was 130 GA/sec. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 23, 1983 
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APPENDIX 


STATISTICAL DATA ANALYSIS 
Probability Plotting Method 

The estimation method used assumes that the parent distribution for 
9 •••/ Y n belongs to some location-scale family. Such a family has the 

representation 


G(y) = G Q ( (y-n )/6 ) 


with density 


g(y) = g Q ( (y-n)/<5) 5 -1 


where G Q (z) is a known distribution, and r\ and 6 are unknown location and scale 
parameters. We refer to the generator distribution G q (z) as a standardized distri- 
bution because it is the distribution of the standardized variable Z = (Y - r])/6 
and serves to specify the location-scale family. 

Let Z be a standardized variable having distribution G Q (z). Let G(y) be 
the true distribution of Y = ln(X). If G Q (z) and G(y) differ only by location 
and scale parameters, then a plot of the theoretical z-quantiles versus the theo- 
retical y-quantiles should yield a linear plot. Departures from linearity indicate 
that G(y) does not belong to the family defined by G Q (z). 

The theoretical versus empirical plot is constructed from the pairs (M^ , Y(i)) 
for i = 1, 2, ..., n, where y ( i ^ y ( 2J ^ ^ y ( n) ^ e orciered sample and 

M-j , , ..., M n are the corresponding V theoretical quantiles of G Q (z). The plotting 

positions suggested in reference 13 are calculated with = G Q ”^(m^), where 


m. = 1 - m (i = 1 ) 

i n 


m = (i - 0 .31 75 ) / ( n + 0.3650) 


(i = 2, 3, .. ., n - 1 ) 


m. 

l 


(0.5) 


1/n 


(i = n) 


The plots shown in figures 2 to 4 are constructed by plotting the ordered trans- 
formed sample values against the theoretical quantiles of the standard normal 
distribution 


G 0 ( z) 


(2tt) -1 / 2 



exp(-w 2 /2) 


dw 
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appendix 


Those plots shown in figures 5 to 7 are constructed with the theoretical quantiles of 
the type I extreme value distribution of largest values, G Q (z) = exp(-e“ z ) for 
< z < <». Whenever a linear plot suggests a good fit to the family of distribu- 
tions defined by G Q (z ), the distribution F(x) = G ([ln(x) - n]/<$) for x > 0 is 
judged to be consistent with the original sample x°, x 2 , . .., x n# For example, a 
standard normal distribution G q (z) leads to a lognormal distribution F(x), whereas 
the type I extreme value distribution G Q (z) gives a type II extreme value distribu- 
tion F(x) • 

The probability plot correlation coefficient (ref. 13) used earlier for testing 
the adequacy for the lognormal model is 


' “ E V(i> / 0 n q 


where 


E”i 


n v 1/2 

2 ' 


and 


1/2 


q = 


-.2 


E <*i - 


With this test the normal distribution (for y 1 , y 2 , y n ) is rejected whenever 

r takes a value in the lower tail of its sampling distribution (r is always posi- 
tive). Percentage points needed to define the upper limits of the rejection regions 
are available in reference 13 for various significance levels. A similar test is 
available (ref. 15) for testing the adequacy of the extreme value model. 


Quantile Estimators 


Let y^ = ln(x^) for i = 1, 2, ..., n denote the transformed sample values 

and let y^ denote the 5th quantile of the true distribution G(y). Asymptotically 
unbiased estimators originally proposed in reference 14 are described below. 

A 

The estimator y^ of y is computed from the ordered sample 
y (1) < y (2) * ••• < y (n) as ^ >the linear function 


y? = c y<M) 


+ (1 - c) 


y (L) 


III III i 111,1 
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where Y(l) < Y( M ) are chosen f rom the ordered sample with indices 
L = [nq] + 1 
M = [np] + 1 


and [np] and [nq] are the integral parts of the real values np and nq. For a 
specified the^quantities 0<q<£<p<1 and 0 < c < 1 are chosen so that 

the variance of is minimized. These quantities are reproduced from reference 14 

in table A1 for the case of a parent normal distribution. 

2 A 

The variance T of the estimator from reference 14 is 

= (6^/n) V(p,q), where 


V( Pr q) = c 2 p 1 ~ ^ + (1 - c) 2 q — ■ *" + 2c( 1 - c)q 

g (z ) g (z ) g (z ) g (z ) 

op o q o p o q 


and g Q (z) is the density function of the standardized variable Z = (Y - n)/6; 
Zp and z^ are theoretical quantiles of the distribution G Q (z). 


In the case of a normal distribution , 


and 


and 


g Q (z) 


is the standard normal density 


are quantiles of the standard normal distribution. For this case. 


~ o 2 

the estimator t of t is computed by replacing 
variance S^, which is computed from y^ , y 2 , • ' 
are shown in table A1 • 


n 


by its estimator, the sample 
Computed values of V(p,q) 


The transformation x^ = exp(y^) converts the estimator y^ of y^ to an 


estimator of the quantile x^ of the distribution F(x). The approximate variance 

A A « A 

of x^ shown in tables 3 to 5 is computed with var(x^) = exp(2y^), with the 

A 

approximate form obtained from a Taylor series expansion of exp(y^). (See method 
discussed on p. 302 of ref. 16.) Confidence limits for x^ follow from the result 

that t ” 1 ln(x^/x^) = t ” 1 (y^ - y^) has a limiting normal distribution (as n ■* ») 
with mean equal to zero and variance equal to one (ref. 14). The upper and lower 
limits U and L for a (1 - a) 1 00-percent confidence interval are 


L = k ex P ( ' z (1-a)/2 x) 
u “ *5 ex P (z (1-a)/2 x) 


where z (^- a )/2 is the u PP er ^ 1 — ot )/2 percentage point of the standard normal 
distribution. 
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Sample Size 

Let R = U/L. From the preceding representation of L and U, ln 2 R is a 
random variable having expected value - a ) ~ 42 ^i - a ) /2 ^ 2 / n ) V(p f q) • The 

sample size needed to ensure that this expected value does not exceed a prescribed 
value ln 2 R Q is the smallest n satisfying 4z ^ (<5 2 /n) V(p,q) < ln 2 R Q . The 
solution is 


n = 4z (1-a)/2 S2 V(p ' q) ln ’ 2R o 

The sample size depends on the index 5 hut not on the parameter x^. It also 
depends on the level of confidence. An estimate of the scale parameter 5 is needed 
to compute n. 

The sample size for a required R increases with increasing £ . This is not 
surprising because the samples yield relatively little information concerning distri- 
butional shape at the tails. For less extreme quantiles (5 < 0.95), the confidence 
limits are narrower and the estimators have good asymptotic relative efficiencies 
(ref. 14). 


TABLE A1 OPTIMAL SELECTION OF INDICES FOR PARENT NORMAL DISTRIBUTION 


5 

c 

(a) 

p 

(a) 

q 

(a) 

z p 

Z q 

V(p,q) 

0.99 

0.302 

0.997 

0.984 

2.75 

2.15 

1 1 .44 

.95 

.313 

.985 

.920 

2.17 

1 .41 

3.64 

.90 

.324 

.970 

.840 

1 .88 

.99 

2 .36 

o 

00 

. 

.371 

.930 

.680 

1 .48 

.47 

1 .63 


a From reference 14. 


9 




1 • 
2 . 


3. 


4. 


5. 


6 . 


7. 


8 . 


9. 

10 . 
11 . 
12 . 

13. 

14. 


REFERENCES 

Thomas, M. E.s Direct Strike Lightning Measurement System. AIAA-81 -251 3 , 

Nov. 1981 . 

Thomas, Robert M. , Jr.: Expanded Interleaved Solid-State Memory for a Wide Band- 

width Transient Waveform Recorder. Lightning Technology, NASA CP-2128, 
FAA-RD-80-30 , 1980, pp. 119-129. 

Von Bokern, G. J.; Piszker, L. D.; and Brick, R. 0.: In-Flight Lightning Data 

Measurement System for Fleet Application. Federal Aviation Administration - 
Georgia Institute of Technology Workshop on Grounding and Lightning Protection, 
FAA-RD-78-83 , May 1978, pp. 345-363. 

Trost, Thomas F.; and Zaepfel, Klaus P.: Broadband Electromagnetic Sensors for 

Aircraft Lightning Research. Lightning Technology, NASA CP-2128, FAA-RD-80-30, 
1980, pp. 131-152. 

Pitts, Felix L.: Electromagnetic Measurement of Lightning Strikes to Aircraft. 

J. Aircr., vol. 19, no. 3, Mar. 1982, pp. 246-250. (Available as 
AIAA-81 -008 3R. ) 

Rudolph, Terence; and Perala, Rodney A.: Interpretation Methodology and Analysis 

of In-Flight Lightning Data. NASA CR-3590, 1982. 

Trost, Thomas F.; and Pitts, Felix L.: Analysis of Electromagnetic Fields on an 

F-106B Aircraft During Lightning Strikes. Proceedings of the International 
Aerospace Conference on Lightning and Static Electricity - Volume 1 , Public 
Relations Section, Culham Lab., Mar. 1982, pp. B3-1 - B3-11. 

Turner, C. D.; and Trost, T. F.: Laboratory Modeling and Analysis of Aircraft- 

Lightning Interactions. NASA Grant NAG-1 -28, Texas Tech Univ., Aug. 1982. 
(Available as NASA CR-1 69455.) 

Pitts, Felix L.? and Thomas, Mitchel E. : 1980 Direct Strike Lightning Data. 

NASA TM-81946, 1981 . 

Pitts, Felix L.y and Thomas, Mitchel E.: 1981 Direct Strike Lightning Data. 

NASA TM-83273 , 1982. 

Thomas, Mitchel E.y and Pitts, Felix L.: 1982 Direct Strike Lightning Data. 

NASA TM-84626 , 1983. 

Wilk, M. B.y and Gnanadesikan, R.: Probability Plotting Methods for the Analysis 

of Data. Biometrika, vol. 55, no. 1, 1968, pp. 1-17. 

Filliben, James J.: The Probability Plot Correlation Coefficient Test for 

Normality. Technometrics, vol. 17, no. 1, Feb. 1975, pp. 111-117. 

Kubat, Peter; and Epstein, Benjamin: Estimation of Quantiles of Location-Scale 

Distributions Based on Two or Three Order Statistics. Technometrics, vol. 22, 
no. 4, Nov. 1980, pp. 575-581 . 


10 



15. Gerlach, B.: A Consistent Cor relation- Type Goodness-of-Fit Test? With Applica- 

tion to the Two-Parameter Weibull Distribution. Math. Operationsforsch. 
Statist., Ser. Statist., vol. 10, no. 3, 1979, pp. 427-452. 

16. Cox, D. R.? and Hinkley, D. V.: Theoretical Statistics. Chapman & Hall, Ltd., 

1974. 


11 


SYMBOLS 


B magnetic flux density, T 

B rate of change of magnetic flux density, T/sec 

B^ rate of change of magnetic flux density for nearby lightning strike, T/sec 

Bg rate of change of magnetic flux density for lightning attachment to 

aircraft, T/sec 

2 

D electric flux density, c/m 

• 2 

D rate of change of electric flux density, A/m 

• 2 

rate of change of electric flux density for nearby lightning strike. A/m 

Dg rate of change of electric flux density for lightning attachment to 

aircraft. A/m 2 


D F' D W,L' D W,R' D T 

b l ,b t ,b w,l ,b w,r 

e x ,e y ,e z 



sign conventions for electromagnetic variables of figure 1 


dc direct current 

E electric field, V/m 

F(x) underlying distribution of X 

G(y) distribution of Y 

G_(z) distribution of the standardized variable Z = (Y - n)/<$ 
o 

g(y) density function of G(y) 

g Q (z) density function of G q (z) 

I current, A 

• 

I rate of change of current. A/sec 

L lower confidence limit 

theoretical quantiles of G Q (y) indexed by integers 
n sample size 

p, q quantile indexes used when referring to the distribution G Q (z) 


R = U/L 

R q prescribed value of R 

r probability plot correlation coefficient 
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u 

V(p,q) 

X 

A 

X C 

Y 


y 

y 


y (i) 

y 5 

A 

y 5 

Z p ,Z q 

1 - a 

6 


n 

5 

2 

T 


variance of the sample y^ , y 2 , • • • , y n 
upper confidence limit 

2 

a factor in the representation of t 

random variable representing a span measurement 

theoretical £th quantile of F(x) 

estimator of 

natural logarithm of X 

value of the random variable Y 

mean of the sample y 1 , y 2 , y R 

ith ordered value in the sample y <1 , y 2 , Y n 

theoretical 5th quantile of G(y) 

estimator of y^ 

theoretical pth and qth quantile of G Q (z) 
level of confidence 

scale parameter of the distribution G(y) 
location parameter of the distribution G(y) 
quantile index 

A 

theoretical variance of y^ 
estimator of t 2 




























TABLE 7.- I FREQUENCY-OF-OCCURRENCE 
DISTRIBUTION 


Lower limit of 
grouping , GA/sec 

Frequency 

2.976 

2 

10.860 

1 

18.750 

1 

26.640 

0 

34.520 

0 

42.410 

0 

50.300 

1 

58.180 

1 

66.070 

1 

73.960 

1 


TABLE 8.- FREQUENCY-OF- OCCURRENCE 
DISTRIBUTION FOR LOGARITHMICALLY 
TRANSFORMED B g SAMPLE 


Lower limit of 
grouping, T/sec 

Frequency 

4.0 

1 

4.4 

1 

4 .8 

0 

5.2 

2 

5.6 

4 

6.0 

6 

6.4 

11 

6.8 

12 

7.2 

9 

7.6 

0 





TABLE 9.- FREQUENCY-OF-OCCURRENCE 
DISTRIBUTION FOR LOGARITHMICALLY 
TRANSFORMED D g SAMPLE 


Lower limit of 
grouping, A/m 2 

Frequency 

0 

0 

.4 

3 

• 8 

6 

1 .2 

1 2 

1 .6 

22 

2.0 

17 

2.4 

24 

2.8 

7 

3.2 

2 

3.6 

0 


TABLE 10.- FREQUENCY-OF-OCCURRENCE 
DISTRIBUTION FOR LOGARITHMICALLY 
TRANSFORMED I SAMPLE 


Lower limit of 
grouping, A 

Frequency 

7.0 

0 

7.3 

0 

7.6 

5 

7.9 

5 

CN 

• 

CO 

8 

8.5 

3 

00 

• 

00 

5 

9.1 

0 

9.4 

1 

9.7 

0 





TABLE 11.- NORMAL PROBABILITY PLOT 
CORRELATION COEFFICIENT r 


Sample 

n 

r 

Significance 

level 

B 

• s 

46 

0.945 

<0.005 

D s 

93 

.933 

>.250 

I 

27 

.980 

>.250 


TABLE 12.- QUANTILE ESTIMATES FOR B g SAMPLE 


5 

A 

x € 

A 

var(x^) 

L 

U 

0.99 

1726.97 

410 369.59 

834.71 

3573.02 

.95 

1726.97 

130 343.63 

1146.40 

2601 .58 

.90 

1497.09 

63 702.34 

1075.83 

2083.31 

• 

00 

o 

1 297 .77 

33 063.77 

986.12 

1707 .90 


TABLE 13.- QUANTILE ESTIMATES FOR D g SAMPLE 


5 

✓V 

x 5 

A 

var (x^ ) 

L 

U 

0.99 

27.33 

36.16 

17.76 

42.07 

.95 

21 .07 

6.83 

16.53 

26.87 

.90 

16.10 

2.59 

13.23 

19.58 

o 

00 

• 

14.06 

1 .37 

11 .95 

16.55 
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Figure 2.- Empirical transformed B g sample versus theoretical quantiles of 

standard normal distribution. 
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Empirical transformed D g sample versus theoretical quantile of 
standard normal distribution. 
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Figure 4.- Empirical transformed I sample versus theoretical quantile of 

standard normal distribution. 
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Figure 5.- Empirical transformed Bg sample versus theoretical quantile of extreme 

value distribution. 
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Figure 6 Empirical transformed D g sample versus theoretical quantile of extreme 

value distribution* 
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Figure 7.- Empirical transformed I sample versus theoretical quantile of 

extreme value distribution. 
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